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Abstract. 

The need for suitable many or infinite fermion correlation functions to describe strongly 
correlated systems is discussed, and the question linked to the need for a correlated 
basis, in which the ground state may be positive definite for certain low dimensional 
geometries. In seeking a positive definite basis a particular trial basis is proposed, 
based on that for hard core bosons in pure one dimensional systems. Single particle 
correlations in this basis are evaluated for the case of the ground state of a quasi- ID 
Hubbard in the limit of extreme correlation. The model is a strip of the 2D square lattice 
wrapped around a cylinder, and is related to a ladder geometry with periodic boundary 
conditions along it's edges. This is done using both a novel mean field theory and exact 
diagonalisation, and the basis is indeed found to be well suited for examining (quasi)- 
order in the model. The model has a paramagnetic region and a Nagaoka ferromagnetic 
region. In the numerical calculation the correlation function in the paramagnetic phase 
has power law decay and the charge motion is qualitatively hard core bosonic. The 
mean field leads to an example of a BCS type model with single particle bosonic long 
range order. 



1. Introduction. 

There is currently a great deal of interest in the question of the movement of charge in 
low dimensional strongly correlated electron systems, partly inspired by the discovery in 
1986 of the high temperature superconductors.^ One of the difficulties in describing such 
systems is the invention and evaluation of suitable correlation functions. The problem 
is to decide what correlation functions most aptly measure the charge movement and 
order (if any) present in the ground state, and the elementary excitations. 

A linked question is, what is the most appropriate basis to use in representing the 
ground state and measured correlation functions? This latter is a complex issue, as the 
most appropriate basis for describing the ground state and elementary excitations may 
not be the most appropriate for describing probes of the system. A clear case in point 
is that of photoemission from a spin-charge separated system. A spin-charge separated 
system is presumably best described using a basis involving separate operators for spin 
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and charge. Photoemission, however, involves the ejection of a complete electron, so 
predicting the photoemission spectrum of a spin-charge separated system would involve 
taking convolutions over separate spin and charge wavefunctions, described using sepa- 
rate spin and charge operators. Clearly not only an appropriate basis is required, but 
also the transformations between it and the bases appropriate to experimental probes. 

In seeking to develop an understanding of low dimensional correlated systems one could 
start by using one or two particle correlation functions, expressed in the most natural 
basis to work in - usually a simple fermionic one. The question is then whether or not 
this is sufficient to describe the behaviour of the system. This is clearly not always the 
case. A simple example is the movement of domain walls in models of magnetism. A 
simple one dimensional Ising model can contain domain walls, which in one dimension 
are single particle excitations of the model. However, creating a domain wall involves 
changing all the spins in the domain - potentially a large number. Moving the domain 
wall from some site i to some other site i+n involves flipping n spins. Hence a description 
in the "natural" basis requires at least n fermion correlation functions. The key to the 
problem is that the motion, viewed in that basis, is collective. However, a collective 
basis also exists, in which we create and annihilate domain walls, using some and b^. 
Expressed in this basis the domain wall motion ia a simple one particle affair, requiring 
only one particle correlation functions. 

The implication is that in describing systems which involve strong correlations and 
hence possible collective motion, we may need to use n or even oo particle correlation 
functions if we wish to use the "natural" basis. Alternatively, we can seek a collective 
basis, in which the correlations become single particle entities. The latter would seem 
the most appropriate choice, at least for describing the ground state charge motion and 
possible order. 

ID systems are known to be spin-charge separated^. In the one dimensional Hubbard 
model, for example, the ground state is uniform phase and the excitations are phase 
excitations, which are indeed collective. In other words, the wavefunction contains 
no nodes, (is positive definite) and excitations can be described as the introduction 
of nodes. The most appropriate basis for describing the ground state charge motion 
in this system is clearly the uniform phase (or "positive definite") basis, otherwise 
known as hard core bosons. The basis is well known from the study of the XY-model. 
The transformation between hard core bose creation and annihilation operators and 
those for fermions (needed for photoemission etc.) is also known - the Jordan- Wigner 
transformation (see section 2). This is a transformation applicable in ID which maps 
spinless fermions (/j, //) onto hard core bosons (6j, bj) and vice versa. In other words, it 
exactly gives the positive definite basis in terms of the fermionic one. It is clear that, for 

example, > (^///j/i+n/j+n) since there will be no phase cancellation 

for the hard core bosons. In other words the positive definite, uniform phase basis is 
also the maximal correlation function basis. This may prove useful in seeking to identify 
the equivalent positive definite basis in quasi- ID and 2D systems. As it is the maximal 
correlation function basis, any ground state order present will be seen most clearly in 
this basis, and order seen in this basis need not be easily detectable in other bases. 

The reason for this interest is that it has been suggested^ that the 2D Hubbard model 
may similarly be spin-charge separated, with a positive definite charge wavefunction. 
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and possibly collective charge motion. How this might come about is still unclear, since 
in a 2D fcrmionic system the nodes are unavoidable. However, one could imagine a 
spin-charge separated system in which all the nodes (the fermi statistics) are swallowed 
up into the spin wavefunction in some way. With a suitable difference in energy scales 
between the two subsystems this could leave an effective node-less charge wavefunction. 
(A similar situation will be examined in this paper.) The question of whether this 
actually occurs still remains contentious. Investigations usually revolve around the 
Hubbard model and it's strong coupling limit, the tJ model. 

Analytic attempts to answer the question have not so far proved conclusive. Some 
renormalisation based calculations^ find fermi liquid behaviour anywhere above pure ID, 
but some other authors disagree.^ Others have attempted to use a numerical approach. 
(See [6] for a detailed review of analytic and numerical work). Results from exact 
diagonalisation (ED) and quantum monte carlo (QMC) simulations yield complex and 
often tantalising phase diagrams, but each have severe limitations. There are indications 
of pairing correlations, usually associated with regions of phase separation at larger J /t. 
Anti-ferromagnetism occurs at half filling, and possibly Nagaoka fcrromagnetism^ near 
to half filling at small J/t^'^~^^. However, the system sizes involved are very small. In 
ED the largest is usually 26 sites^'^, while QMC simulations^'^" can manage slightly 
larger systems, arounnd (12x12) but results are dogged by the sign problem^ and cannot 
reach low temperatures. Hence, direct ED and QMC calculations on the 2D Hubbard 
and t J models are still too restricted, and finite size effects make the boundary conditions 
far too important to draw solid conclusions. In particular, there are no clear indications 
that the pairing correlations lead to superconductivity in either the Hubbard or the tJ 
models.^ Nor is it clear whether the excitation spectrum is essentially fermionic or spin 
charge separated. More promising, perhaps, is the recent 2D momentum space DMRG 
calculation by Xiang^^ which almost exactly matches the accuracy of ED calculations at 
smaller system sizes (< 26 sites) and can also produce results on system sizes of (to date) 
(12x12). Quantitatively, results are in keeping with QMC, but are now true variational 
bounds without sign problems. Further work along this line may be promising. 

Another line which also seems promising is the study of doped tJ and Hubbard models 
on ladder geometries - see review in [14]. Whilst these have direct relevance to ladder 
materials, such as SrCu20^ and Sr2Cu.^O^, they can also be viewed as finite width 
strips taken from the 2D square lattice. The finite width means that charges can move 
around one another, so exchange statistics start to have some importance, unlike in 
pure ID. Hence one might seek to use these models to infer the behaviour of the 2D 
models, by extrapolating results from either single strips^^'^^ or from a series of wider 
and wider strips. Reasonable results are availible for (at least) 2 leg,^^~^^ 3 leg-'^^'^^ 
and 4 leg^^'^^ ladders. As one would expect from renormalisation group arguments, they 
all spin-charge separate, but the excitation spectrums vary massively : systems with 
odd numbers of legs form Luttinger liquids, with gapless excitation branches for both 
spin and charge. Even leg systems, on the other hand, are gapless only for charge 
excitations, and have been characturised^"^'^^'^"^ as Luther-Emery liquids. This large 
difference makes extrapolating the 2D behaviour almost impossible. The problem can, 
in a sense, be blamed on boundary conditions and finite size effects again : the edges of 
the strip constitute an open boundary condition running the length of what is in reality 
a very narrow system. 
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This last point is one that we will address in the current work. The model we will study 

(q.v.) can be viewed as a ladder geometry with periodic boundary conditions along the 
edges. As such it is hoped that it will prove more useful than the ladders in predicting 
the 2D behaviour. 

In seeking an answer to the question of spin-charge separation in 2D one must also 
find a suitable basis in which to describe a 2D spin-charge separated system with a 
positive definite charge wavefunction. Currently, we do not know what this basis is 
(though, some mathematical progress will be made in a subsequent paper^^,) and so 
we will NOT provide any kind of an answer to that question. Instead, we will simply 
move towards a 2D description, by looking at a particular strongly correlated quasi-lD 
Hubbard model, which does spin-charge separate, introduced in section 3. First, in 
section 2, we introduce the basis and correlation functions we intend to use. We wish 
to seek a positive definite basis for our quasi-lD system, which we anticipate will be 
collective, and to examine single particle correlations defined in this basis. Unfortunately 
we do not know what such a basis should be, so instead we use a spin-charge separated 
basis introduced earlier^^. We then apply the Jordan- Wigner transformation to the 
charge subsystem, and use the resulting basis as a start point. The ground state charge 
motion and order are examined in this basis using an exotic mean-field technique in 
section 4, and exact diagonalisation in section 5. In section 6 we conclude. 



2. The Basis And Correlation Functions. 

Before we can write down the correlation functions, we must define the basis. The 
normal start point for fermionic systems is the fcrmionic basis //^, ^. We wish to 
restrict ourselves to the limit of infinite Hubbard repulsion (U — > oo) so we project 
out the doubly occupied states :- 

gla = fl (i - fUi,,) 

The spin-charge separated operators^^ may be written as 

where c\_^ and c^^ are creation/annihilation operators for a single spinless fermionic 
charge at site i^. The phase is chosen to be positive for wavefunctions in which the 
fermionic operators occur in the order Cq, c|, c^, cj^, ... The spin degrees of freedom 
are described by a compressed spin chain \a]^a2----cr^ acted on by ring exchange oper- 
ators R^^ Q,^. Note that ^ Z labels the position of spin = ±^ in the compressed 
spin chain, not the actual lattice site on which it occurs. -R^^^ cyclically permutes all 
the spins from to a^. The ring exchange operators are related to the more common 
Heisenberg spin operators S by the relation 



f3=ai 



(3) 
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Any electron wavefunction can of course be expressed in terms of these operators, be 
it fermi liquid, or spin-charge separated (Luttinger or Luthcr-Emcry liquid) or other- 
wise. However the operators are clearly more suited to describing spin-charge separated 
situations, and are most easily applied in ID or quasi-lD where an obvious linear site 
labeling exists. An exact mapping between these new operators and the conventional 
operators was given in reference [25]. The operation of this mapping can be shown by 
example. The state 

^ ^ ^ ^ written as 9l,^9l,i9l,^\0) (4a) 
in the standard notation, is now expressed as 



..©... 



written as CIC3C4I Tii2T3)- (4b) 
Similarly the operator 

Yl9la9i+2,a becomes J24c^+2[{^ - cl+i^i+i) + 4+lCi+lRa+l,a\ (5) 

i [i:a] 

The first term in square brackets deals with the situation in which there is no fermion at 
{i + 1) and hence the order along the compressed spin chain is unchanged. The second 
term deals with the situation in which (^ + 1) is occupied so that two spins in the chain 
are cyclically permuted. The summations over i (lattice sites) and a (positions in the 
compressed spin chain) are not, of course, independent and so care has to be taken in 
performing them. 

This spin-charge separated basis is the one in which the model is most easily written 
(q.v.) and in which the mean field and numerical solutions will be most easily expressed. 
It is not, however, the positive definite basis we seek, in which we wish to describe the 
ground state charge motion. The charge degrees of freedom are here still expressed in 
terms of spinless fermion operators. Exactly what the required basis is for this model 
we do not know. As an approximation to it we will use the basis given by the Jordan- 
Wigner transformation.^^ This can be expressed as 

^'i = c^exp ^iTT^cjc^.j (6) 

where 6| and create and annihilate hard core bosons on the site z, and the phase 
is again chosen to be positive for wavefunctions created using fermionic operators in 

the order Cq, c|, C2, .... , cj,^, The itt phases from the exponents account for the 

difference in exchange statistics between fermions and bosons. Hence, for example, 

= 4+2Ci (1 - M+l^i+l) (7) 
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This gives us a spin-charge separated basis in which we hope we have something close 
to a uniform phase representation for the charge wavefunction. It is approximation, 
however, because the Jordan- Wigner transformation works exactly only in srict ID. 
The charge correlations we wish to measure are then single particle correlations in this 
basis - 

Bn = {blb,+n) (8) 

These may turn out to be not quite maximal, as this is not the pure ID chain for which 
Jordan- Wigner is defined. 

In terms of the original basis they are also collective. This can be seen by back trans- 
forming - 




Hence, written in the spinless fermionic basis this is a complex n-particle correlation 
function, describing the movement of a charge over n sites, subject to the positions of 
the charges in between. However, in the Jordan- Wigner basis it measures the mobility 
of what, if this was a pure ID system, would be hard core bose quasi-particles. Al- 
ternatively it can be viewed as measuring the movement of a bare charge over n sites, 
irrespective of fermi exchange statistics. If is long ranged 7^ 0) in a partic- 

ular ground state it means that thereare long ranged charge fluctuations, but this does 
not ncccasarily mean that the state is superconducting. This is because, flrstly, we will 
present no information on the excitation spectrum of the model, so we can only hope 
that the low lying charge excitations are gapless, with the same statistics and mobility. 
Secondly, even if that is the case, we also require a macroscopic number of charges to 
be involved. However, if j!^!^B^ is zero it does mean that there is no order, and that 
the state is not superconducting. 

If the correlation function in (8) is applied to a non-interacting free electron gas it is 
found to decay^^ with the power law ~ -\. 

n2 



3. The - tg Model. 

We are interested in the 2D square lattice, in the extreme strongly coupled limit, U = 00, 
J = 0. As mentioned above, a popular approach to the problem of the 2D square lattice 
is the study of (flat) ladder geometries. Some alternatives were introduced in reference 
[28], having different boundary conditions, and so providing alternative scalings to 2D. 
The one we study here is equivalent to wrapping a section of the 2D square lattice (with 
nearest neighbour bonds only) around a cylinder. It can be viewed as a multi-leg ladder 
with periodic (instead of open) boundary conditions parallel to the main axis. We also 
introduce a one lattice site displacement at this boundary. This results in a single spiral 
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of sites running along the cylinder, (see figure 1,) with two distinct types of bonds, 
one around and one along. (These may be viewed as to some degree analagous to the 
perpendicular and parallel bond directions in the ladder geometries.) 

We elect to label our sites along this spiral. This gives us a useful ID labeling scheme for 
our quasi-lD model, in which each site is connected to two nearest neighbours, and to 
two n^^ nearest neighbours. It should be noted, however, that both of these arise from 
nearest neighbour bonds of the 2D square lattice {not next nearest neighbour bonds) 
and so, ideally, we would like to solve the models with the hopping parameters equal in 
these two directions. 

What we have described is a homologous series of t^-t^ models, with n odd. n is 
proportional to the diameter of the cylinder, so the n — > oo limit recovers the 2D 
lattice. Note that it is the bipartite nature of the 2D square lattice that forces n 
to be odd. Consideration of longer range hopping in pure ID leads naturally to the 
consideration of the ^1-^2 model. However, we are interested in 2D, not ID, and t^- 
t2 actually corresponds to a strip from the 2D triangular lattice, where non-bipartite 
frustration leads to very different physics from that in which we are interested. Thus, 
the first model in the series is t^-t^ , involving first and third nearest neighbour hopping 
terms at each site :- 

H = -h J2^9la9r+l,a + ^c] - [^7j,,^7,+3,a + C.C.] (10) 

Since all the t^-t^ models are one dimensional in a renormalisation group sense they 
will all exactly spin-charge separate. Written in terms of our spin-charge separated 
operators t-^^-t^ becomes 



^q+i + c.c. 



tl,t3 



i 

-h Yl [4cj+3(l - 4+iq+i)(l - c|+2Ci+2) + C-C. 

i 

[4*^1+3 (^^i+l^i+li^ ~ ^1+2^1+2) + (1 ~ ^i+l^i+l)^i+2^i+2^ ^a+l,a + 

[i:a] 

\_^i^i+3^i+l^i+l^i+2^i+2-^a+2,a + 

[i-.a] 

(11) 

In terms of the spiral labeling the first term is simply nearest neighbour hopping and 
so involves no change in the compressed spin chain. The other terms are for third 
neighbour hopping over zero, one and two intermediate charges respectively. Hence the 
second term involves no ring exchange, whist the third and fourth terms involve ring 
exchanges of two and three spins respectively. 

The t^-t^ series of models has been examined previously^^'^^ when an exact perterbation 
theory in the limit ^ — was presented. This calculation was for an infinite chain at 
T = OK. For the t^-t^ model a paramagnetic state was found below a critical fiUing^^ 
n^ = 0.6675045, and a Nagaoka^ ferromagnetic state above it, in keeping with some 
of the numerical results for the 2D square lattice^, for example. It should be noted. 
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however that in this instance (unhke in the numerical results) this is the Nagaoka 
state in the thermodynamic limit of a finte hole density. It was found, however^^, 
that the value of n^ was dependant upon n, and that n^ ^ 1 as n — > oo. le the 
ferromagnetic phase disappeared in the limit of 2D square lattice connectivity, (albeit 
with the "wrong" matrix elements.) Unfortunately we are unable to prove this away 
from the ^ — > limit. Here we will present evidence that the phase diagram of t^-t^ 
at least is qualitatively unaltered away from this limit, and that the same phases are 
found even in the physically interesting case of |^ =1. 



4. Mean Field Treatment Of The - Model. 

The basic principle is to study the spin degrees of freedom in the mean field of the 
charge degrees of freedom, and the charge degrees of freedom in the mean field of the 
spins. This amounts to assuming that the spin and charge are not only separated (as 
expected in ID) but "independent." In this way local correlations between the two are 
lost, in favour of an effective spin hamiltonian and an effective charge hamiltonian. In 
the limit 1^ — ^ ^^he mean field recovers the exact perterbation theory results^^. The 
charge hamiltonian was 

i 
i 
i 

"-^2^3 Y2i [4<^i+34+l<^i+l4+2<^i+2 + ^-cj 



charge 



C.C. 



where 



The spin Hamiltonian was 



(12) 



Kin = [< + ^^2 + ^3°] K + [N^ + ^3 ] E + ^-c] + iV| J] [K+2,a + CC 

a a 

(14) 

is the expectation value of the operator which moves a charge from i to i + n passing 
m occupied sites on the way. To express it in terms of the c| operators, we first define 
the occupancy of a site j as Xj = or 1. Then the set {Xi+i,Xi+2^ -^Xi+j, •••,Xj+n-i} 
indicates which of the sites between i and i + n are occupied. Hence Yl^=i Xi+j — 
The operator which checks whether or not a particular site i + j is occupied in the 

wavefunction is then ^1 — c|^.jC^^.j-^ (cl_^_jC^_^j^ . Finally, in evaluating A^^^ we 

need to sum over all possible permutations of the set {..,x^-,..}, ie we sum over all 
possible ways to have m sites occupied between i and i + n. Writing the sum over 
permutations as } we obtain 
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n 



^n{^l^i+n^i,n) 



(15) 



where 



i,n 



n-1 



P{{Xl,X2,---,Xj,---,Xn}) J = l 



n-1 

(16) 



The two sets of equations were examined self consistently as a function of the filling 
TT-g and the ration ^ . The spin hamiltonian is actually the well known J1/J2 nearest 
and next nearest neighbour Heisenberg model. The charge hamiltonian is a little more 
complex, and requires a further approximation. A Bogoliubov-Hartree-Fock mean field 

approximation was used, in which pairing correlations, ^c|c|_|_^^, were allowed. This 

led to a BCS type model with charge pairing, which was solved using a Bogoliubov 
transform^^. A variational phase diagram was found by comparing the energy of the 
charges subject to various trial spin wavefunctions. The spin states used were the Neel 
state, the Heisenberg ground state, and a ferromagnet. For the low spin states the 
ground state always involved pairing. Initial results were published in reference [30]. 
The phase diagram obtained is shown here in figure 2. Qualitatively almost nothing 
alters as we move away from the |^ — > limit. It was previously noted that the drift 
of the phase transition towards lower fillings as |^ increases is probably due, at least 
in part, to the fact that the spin-charge separation mean field is variational for the low 
spin states, but exact for the ferromagnet. 



Here, in an extension to the previous paper, the results are presented for the correlation 
function introduced in equations (8), (9). Since this actually measures correlations 
in the ID hard core bose (Jordan- Wigner transform) basis it is worth noting that the 
effective charge hamiltonian becomes the hamiltonian for hard core bosons (XY model) 
on the t^-t^ lattice if we set = —1 and R2 = +1. For these values, measures 
maximal correlation functions by definition. This will be used for comparison shortly. 



In the mean field. 



B^ becomes^^'^^ 



n-1 



a. 



a 
a-, 



1 

2 



So 



a. 



2 

n-1 



ai — 61 



«n-2 + ^r, 
«n-3 + ^r. 
«n-4 + 



«o - 2 
ai — 5i 



e 

V 

e 
n 



(17) 



Here 



and 



(18a) 
(186) 
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Both can be calculated easily once the gap equations have been solved self consistently 

and the charge ground state found. The |-.|g„g^ indicates that only terms containing 
even powers of 5's (such as 6f or 6f62) are included in the determinant. This can easily 
be achieved by adding two determinants together in which the signs of the 5's are all 
reversed. Hence 



5n = 2 



n-2 



+ 



+ 5, 



«n-2 
«n-3 + ^ 



n-2 
n-3 



«i + ^1 "o - i 



Q!^_2 - 



Figure 3 shows evaluated for the hard core bosons, as a function of n, for ^ = 1 
and filling Ug = 0.5 (quarter filling). The correlations tend to a finite limit of 
= 0.2284785 ± 0.0000001 as n — oo. Also shown on the same axes are the more 
conventional correlation functions = ^cjc|^i^ = 0.1889107 ± 0.0000001, and C„ =< 

cJc|_|_]^Cj_)_^Cj_|_^_j_^ >. The latter is given in mean field by 

C^ = 5f - al + q;„_iQ!„+i (20) 

and tends to the value of as n — > oo. Being maximal for hard core bosons, B^ is 
indeed the largest. 

Figure 4 shows J^l!^B^ versus n^ at |^ = 1 for hard core bosons and two paramagnetic 
variational states:- charges moving subject to the Heisenberg and Neel states along the 
compressed spin chain. It should be noted that for charges moving in the ferromagnet 
there is no scaling of 5^, which occillates wildly in sign. This is not suprising as the 
charges here are known to be spinless fermions. For charges moving subject to the 
Heisenberg ground state along the spin chain (the variational ground state for |^ = 1.0 
and Ug < 0.587) there is again no scaling behaviour for B^ above ng=0.5939, and the 
gap [5]) closes very rapidly at this point. This behaviour will be easier to explain in 
the context of the numerical calculations, q.v. It is not a problem here, however, as 
it occurs just above the point at which this state ceases to be the variational ground 
state. (Unfortunately problems with the numerics make producing accurate data for 
0.5639 < Ug < 0.5939 impossible.) 

In figure 5 the values of J^^l^B^, n^-^Ioo^n plotted together for the variational 

ground state (Heisenberg) at ^ =1.0. For the Neel case they are similar. Non-zero 
values (Vn) of B^, and all pairing correlations 5^, are seen for both low spin states, 
indicating that, for most values of Ug, charge pairs exist in the ground state, and have 
long range order. The order is of course a result of the Bogoliubov-Hartree-Fock mean 
field making it an independent charge model. Hence, a "large" number of charges can 
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move freely through the lattice in the mean field ground state. This does not make the 

model a superconductor as the model is one dimensional, and hence the Mermin- Wagner 
theorum forbids a macroscopic number of charges being involved. It should be noted 
though, that, potentially, only the dimensionality prevents macroscopic numbers being 
involved. 

As anticipated, is always the largest of the correlation functions measured. This 

suggests either that 6^ is indeed the maximal correlation function basis, (and hence 
the positive definite one,) or that it is at least closer to it than the conventional basis 
involving Cooper pairs c|cj_|_^. 

Finally, it should also be noted that this demonstrates that long ranged single particle 
charge fluctuations can occur in an BCS type model. This raises the question of whether 
or not a single particle description (in a suitable basis) can be found for other BCS states, 
and if so, would it make a better description of the ground state than the standard 
ones^^? 



5. Exact Diagonalisation Studies Of The - Model. 

The ^1-^3 model has also been examined using exact diagonalisation of 12 and 16 site 
chains, for a variety of electron fillings, as a function of |^ . Periodic or anti-periodic 
boundary conditions were used, as appropriate to the electron filling. The exact di- 
agonalisation was done using the Lanczos algorithm, on SUN SPARK 1, 2 and 10 
workstations. The first results obtained were for the total spin of the system at each 
set of values considered. The results are given in figures 6a and 6b. The total spin 
corresponding to ferromagnetism {S^^^) varies with the number of electrons involved, 
(N^) and if there is an odd number of electrons then there is always a spinon left over 
even in the lowest possible spin configurations. This is awkward to plot, so the quantity 
actually plotted is 

( = for even 

^scaled ~ \ o i „„ i (^"'"^ 

^ gmax_l JVe — 1 6 

This gives a value ranging from for a pure paramagnet to 1 for a ferromagnet. 

The structure of these results is very similar to that of figure 2 (the mean field results) 
with a ferromagnetic region near to half filling, and a paramagnetic ground state at 
lower fillings. The main difference is that here the transition does not sink to lower 
fillings as |^ grows. (The non-zero spin state around = 0.5, |^ = 1.0 is perhaps a 
finite size effect, or possibly a commensur ability effect.) This indicates that the increase 
in the size of the ferromagnetic region seen previously was indeed purely a results of the 
approximation. 

Since the results for 12 and 16 site systems are qualitatively identical, and quantitatively 
almost identical, only the results for 16 sites will actually be shown from this point. 

Following the spin state results, charge-only wavefunctions were extracted from the full 
electron wave-function obtained in the exact diagonalisation. For a finite sized system 
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this can be done whether the system is spin-charge separated or not. Most importantly, 
it is possible to do it in such a way as to exactly obtain the positive definite basis for 
this charge system. The amplitude A^^,^ used for the n charge configuration {Xi} was 
the root sum of squares of the amplitudes 0'{xi},{aj} each of the states in the full 
fermionic wavefunction having that particular charge configuration. le, 

i ai a j an i 

{g\ as defined in equation (1).) The choice of phase is arbitrary, but always choosing 
positive phase gives the positive definite basis. It is possible to do this exactly as we have 
only a finite sized system. However, limitations on system size mean that at this point 
we can extract little useful scaling information for directly. Instead, we first find 
effective charge models and then scale those. This is done by overlapping the extracted 
charge wavefunctions with the ground state of an effective charge model, and varying 
the parameters of the effective model to obtain the optimal fit. The effective model 
used was the charge model that arose in the mean field calculation 

-^i + C.C. 

i 
i 

~^1^3 Z [c|c^_,_3 ^c|_|_;^Cj_|_]^(l — c|_|_2Cj_|_2) + (1 — cJ_|_iC^_|_i )cJ_,_2Cj_|_2 j + C.C. 
[i:a] 

~^2^3 Z [c|c^_|_3c|_,_]^Cj_|_]^c|_|_2Cj_|_2 + C.C. 
[i:a] 

(23) 

Appropriate choice of the parameters can make this an interacting or non-interacting 
spinless fermion model, an interacting or non-interacting hard core boson model, (XY 
model,) or anything in between. This includes charges moving subject to a Neel an- 
tiferromagnet or to a Heisenberg ground state along the compressed spin chain. The 
effective model ground states were again found using the Lanczos algorithm. The size 
of the optimal overlap is indicated in figure 7. The overlaps are generally very large, 
but do become smaller - to at worst 97.3% - near the phase transition itself. 

The effective models can be classified by their proximity to one of the four exactly known 
points, as shown in figure 8. The four points, labeled A, B, C and D in the figure are :- 



A 


Spinless fermions. 




= +1, i?2 


= +1, 


B 


Hard core bosons. 


Ri 


= -1, i?2 


= +1, 


C 


Charges in Neel state. 


Ri 


= 0, i?2 


= 0, 


D 


Charges in Heisenberg 


Ri 


= 1 - 21n2, i?2 


= l-61n2 + |C(3) 




ground state. 




= -0.386294361 


= -0.454255051 



where C,{n) is the Ricmann zeta function. (The values for the Heisenberg case come 
from [32] and [33].) Anything away from points A and B can be interpreted either as 
interacting spinless fermions or as interacting hard core bosons. 

Using this classification scheme a phase diagram can be built up for the behaviour of 
the charges in the ground state of both 12 and 16 site systems, as shown in figure 9. As 
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an example, a few points are indicated on figure 9 and their positions on axes R2 vs 
are shown on figure 10, which should be read together with figures 8 and 9. Note that 
those classified as Heisenberg like are in fact often much more efficient than the charges 
moving in the Heisenberg ground state itself. By increasing the strength of the short 
ranged correlations Ri and i?2 from around 0.386 and 0.243 to ~ 1.0, (presumably at the 
expense of higher R^ correlations not involved in the Hamiltonian) they achieve better 
energies. It should also be noted that the charge motion in the entire paramagnetic 
region is qualitatively that of hard core bosons. 

Once the effective models had been obtained they were themselves exactly diagonalised. 
This could be done for a variety of larger system sizes as there are fewer degrees of 
freedom in the effective charge-only models. This meant that finite size scaling of the 
results could also be done. The quantity evaluated was again the correlation function 
introduced in section 2: 

Bn = (blK+n) (24) 

This was done using the Jordan- Wigner transformation again. The effects of the periodic 
boundary conditions on this number complicate matters. For an infinite ID system, 
is expected to have power law behaviour. For a finite sized system with periodic 
boundary conditions the correlations never completely decay, so the system is ordered, 
and the order propagates both ways round the ring. We expect something like 

5^ = ^ (n-" + (M - n) (25) 

where M is the number of sites on the ring, A is some constant and a some power. 
We could now fit a{M) using (25), then scale with M, but this would be error prone. 
Instead we note that selecting just n — ^- (M even) reduces (25) to 

5n = 2A(^-j (26) 

This can be scaled with M fairly simply. A plot of (S^n*^""^-*) against ^ is used. For 
the correct value of o: a line of zero gradient is obtained, which will show up at larger 
n, irrespective of finite size effects at smaller n. Away from this value of a the curve 
diverges or goes to zero. A variety of trial values are used - see for example figure 11, 
for pure hard core bosons at quarter filling. A value of a = 0.225 ± 0.025 is obtained. 

Table 1 shows the results of a for the linear chain {t^ = 0), and (with ^3 = 1.0) hard 
core bosons and charges in the Neel or Heisenberg ground states. If is calculated 
directly for spinless fermions the sign fluctuates wildly, with no scaling, just as in the 
mean field. This again indicates that Jordan- Wigner does not give a positive definite 
wavefunction for spinless fermions. This does not mean such a basis doesn't exist of 
course, indeed for a finite sized system it must. If a suitable basis could be found, power 
law decay of a single particle correlation function in it could perhaps be expected, albeit 
with a large power. 

Similarly, there is again no scaling behaviour for charges moving in a Heisenberg ground 
state for ng>0.5 and for |^ ~ 1.0. For the Heisenberg ground state, R^ has the same 
sign as for hard core bosons, but i?2 has the opposite sign. The point at which the 
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scaling behaviour is expected to cease is when the Jordan- Wigner basis used for 
ceases to be quahtatively similar to the positive definite basis sought, adn hence the 
wavefunction has positive and negative parts of the same order. This in turn is the point 
at which hopping over 2 other charges, with the phase of i?2, becomes more important 
than hopping over 1, with the phase of R^. This should indeed occur at rig = 0.5. 

Figure 12 gives the power laws obtained from scaling the effective models themselves. 
Points for which effective models were scaled are marked by a cross. Power law scaling 
of is found everywhere except where there is ferromagnetism and hence spinless 
fermionic charges. For all other values of the parameters power laws are found, empha- 
sising the essential unity of the whole paramagnetic phase. The existance of power laws 
suggests that power law behaviour for the charge correlations might also be expected 
in the full quasi-infinite system. This is in keeping with the Hartree-Fock mean field 
results in section 4, except that the spurious long ranged order does not occur in the 
exact diagonalisations. 

In some areas of the |^ versus n^ phase diagram (figure 12) the exponent in the power 
law for the effective model is actually smaller than that expected for hard core bosons. 
This means that the charge correlations found in the ground state of the effective model, 
and hence most likely in the original fermionic model, are longer ranged than in the t^-t^ 
hard core boson model. The longest range correlations occur at the same place that the 
overlap is smallest. This is to be expected, since the effective models used involve only 
local charge interactions, and will therefore have the most trouble describing situations 
where the range of the charge correlations is greatest. 



6. Conclusions. 

In seeking a better understanding of uniform phase bases for strongly correlated sys- 
tems we have introduced a certain spin-charge separated basis. Since we do not know 
how to represent positive definite wavefunctions in (l+S) D or above wc have tried to 
represent the charge wavefunction using operators arrived at using the Jordan- Wigner 
transformation. This would produce the desired positive definite representation if the 
system were in pure ID. We have evaluated the single particle correlations B^ in this 
basis - collective correlation functions in the "natural" fermionic basis. We have done 
this for a particular quasi-lD Hubbard model, the ^1-^3 model, using an exotic mean 
field theory, and exact diagonalisation studies. 

The mean field theory produces a BCS type model, with two ground state phases. Close 
to half filling a Nagaoka ferromagnetic state was found. At lower fillings a paramag- 
netic state was found, in which, in the mean field, there were gaps to single particle 
fermionic excitations, and Jj^-B^ was finite, indicating the existance of order. It was 

also larger in magnitude than the (^clcl_^_^^ type pairing correlation functions consid- 
ered conventionally. The exact diagonalisation studies produce the same picture, but 
without the spurious long range order permitted by the mean field. The charge motion 
in the paramagnetic region was seen to be qualitaively that of hard core bosons, and 
B^ was found to decay with (measured) power laws. The phase diagram is similar to 
those proposed for the 2D square lattice Hubbard and tJ model in refs [6,8]. 
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The ^1-^3 model thus appears to form a quahtatively hard core bose charge ground state, 
with long range collective charge fluctuations. This does not constitute a macroscopic 
superconductor, as the Mermin- Wagner theorum does not allow order in ID, and we also 
know nothing from the current work about the low lying modes. As a result we cannot 
say whether to classify the t^-t^ model as a Luttinger or a Luther-Emery liquid. However 
we can say that the uniform phase charge wavefunction seems to arise as a result of 
the spin-charge separation : the nodes arc placed in the spin wavefunction, resulting in 
strong short range spin correlations which effectively absorb the fermi statistics. The 
numerical results indicate that this is done extremely efficiently. 

This basis is probably not actually maximal, (equivalent to actually being the positive 
definite/uniform phase basis). However, since the correlations measured in it are the 
largest we can find we conclude that it is at least the closest to maximal that we have 
availible. In the case of the numerical results we do come closer - we have the positive 
definite basis exactly - but we cannot evaluate the order in this basis directly, as we do 
not have a single (or few) particle formulation of it. 

The t^-t^ model is a concrete example of a model in which strongly correlated charge 
motion requires many fermion correlation functions for it's description. Here we man- 
aged with an n particle correlation function. In the general case an infinite particle 
correlation function might be required. However, we have also found that single parti- 
cle correlations in a suitable positive definite basis provide a better description. This 
clearly has implications for proposed 2D spin-charge separated systems. 

The ^1-^3 charge dynamics also form a concrete example of a BCS type model for 
which single particle correlation functions can be used in describing the ground state, 
as opposed to the two particle description normally used. This raises the question 
of whether or not a single particle order parameter could be used to describe other 
superconducting states. Half an answer is given - if 5^ is non-zero there may be a 
superconducting ground state (subject to the involvement of macroscopic numbers etc.) 
On the other hand, if — there will not be order and superconductivity. 

None of our results tell us anything about the groundstate or excitation spectrum of 
the 2D square lattice Hubbard, tJ or even t models. However, it is notable that the 
mean field and exact diagonalisation phase diagrams for ^^-^3 away from ^ = are 
qualitatively and almost quantitatively the same as the exact perturbation theory results 
at 1^ ^ 0. In the n — > 00 limit of the perterbation theory, (corresponding to 2D 
connectivity,) only the paramagnetic phase remained. One could suggest from this that 
the 2D square lattice t model should have nothing but paramagnetism anywhere away 
from half filling, with qualitatively hard core bosonic charge correlations, and maybe 
ordering. This would certainly be in keeping with much of the other work on the square 
lattice and ladder geometries, but this is pure conjecture. To proceed beyond conjecture 
we must apply these ideas above ID. The problem with this is one of formulation. Here, 
power law order in ID was seen (suggesting the possiblity of full long range order in >2D 
or perhaps >2D) because something close to the transformation to a positive definite 
basis is known (the Jordan- Wigner transformation). If suitable transformations could 
be found for other systems then perhaps full single particle long range order could again 
be seen. This would obviously facilitate the development of a new way to visualise and 
describe superconductivity in general, and would have wide-ranging application to the 
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study of high temperature superconductivity in particular. 
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Table 1. 

Table showing powers a for scaling of for the specified charge models. Results given 
±0.025. 





a For Power Laws At Given Fillings. 
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Figure 1. 

The first three ti-t^ models, (n = 3, 5, 7) 
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Figure 2 

1^ vs. rig phase diagram for the t-^^-t^ model from the mean field theory. From ref 
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Figure 3 

Scaling correlation functions with n, for hard core bosons. Filling = 0.5 and |^ 
A: = (blb.^J), B: S, = (clcl,), C: Sf, D: < 44+iC,+nCi+n+i >■ 
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Figure 4 



-Jim—B^ as a function of filling n^ at |^ = 1. 
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A: Hard Core Bosons and charges in B: the Neel state and C: the Heisenberg ground 
state. 
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Figure 5 

Correlation functions vs. Ug for mean field variational ground state. (Charges in the 
Heisenberg ground state.) 
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Figure 6a 

'^scaled ^P^^ phase diagram from exact diagonalisation of the t-^-t^ model on 12 sites. 
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Figure 6b 

'^scaled ^P^^ phase diagram from exact diagonalisation of the ti-t^ model on 16 sites. 
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Figure 7 

Error in effective model: (1 - the overlap) between the charges in the full fermionic 
ground state and in the effective charge model, on 16 sites. 
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Figure 8 

Classifications used for the effective charge models with and i?2- 
"x" refers to an exactly known point - see main text. 
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Figure 9 

Phases in the fermionic model as characterised by the effective charge model. 
The selected sample of points corresponds to those shown in figure 10, where their ef- 
fective Ri and i?2 values are plotted. 
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Figure 10 

Effective charge model classifications of the points selected in figure 9. The points 
marked by circles lie along |^ =1.0 in figure 9, and those with crosses lie along |^ =0.1. 
The Ug values are labeled. 
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Figure 11 

Scaling to find the power a. for pure hard core bosons at ng = 0.5, for the tyt^, model. 
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Figure 12 

Powers a for the decay of 5^, from scaling the effective model. Plotted against ng and 
the ratio |^ for the original model. Results are ±0.025. 

The crosses mark the points at which a was evaluated. Note that the exact results^ 
gave a = 0.5 everywhere for |^ =0.0. Near Ug = 1.0 the ground state of the model 
becomes fermionic and hence B„ does not scale. 
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